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Abstract — The optimal PMU locations to collect voltage phase 
angle measurements for detecting line outages in wide-area 
transmission networks are investigated. The problem is estab- 
lished as one of maximizing the minimum distance among the 
voltage phase angle signatures of the outages, which can be 
equivalently formulated as an integer programming problem. 
Based on a greedy heuristic and a linear programming relaxation, 
a branch and bound algorithm is proposed to find the globally 
optimal PMU locations. Using this algorithm, the optimal trade- 
off between the number of PMUs and the outage detection 
performance is characterized for IEEE 14, 24 and 30 bus systems. 
The algorithm is shown to find the globally optimal PMU 
locations in a small number of iterations. It is observed that 
it is sufficient to have roughly one third of the buses providing 
PMU measurements in order to achieve the same outage detection 
performance as with all the buses providing PMU measurements. 

Index Terms — Phasor measurement unit, location selection, 
outage detection, transmission networks, branch and bound 



I. Introduction 

In high voltage transmission networks, lack of wide-area 
situational awareness has been one of the major causes of 
large-scale blackouts Q). One of the reasons for the lack of sit- 
uational awareness has been the limitations of the conventional 
sensors and SCADA (supervisory control and data acquisition) 
systems. Phasor measurement units (PMUs), as compared to 
conventional sensors, are able to provide GPS-synchronized, 
more accurate and temporally much denser measurements of 
voltage and current phasors 0. Thus, the deployment of 
PMUs for grid sensing is widely considered to be a major 
driving force for improving the reliability of transmission 
networks. 

There has been considerable research investigating how 
PMU measurements can be exploited in various tasks for 
achieving a reliable grid, including outage detection O, H, 
0, state estimation 0, Q, O, stability analysis 0, etc. 
As the cost of installing and networking PMUs is relatively 
high, two of the major questions are i) what is the necessary 
or sufficient number of PMUs to use for achieving good 
performance in all these tasks, and ii) where are the optimal 
locations to collect PMU measurements? For achieving full 
network observability, it was estimated that about one third of 
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the buses need to provide PMU measurements JTJ. However, it 
is unclear whether this estimate applies to other kinds of tasks 
including outage detection. Moreover, for different tasks, the 
objectives for which the PMU measurements are used differ 
considerably, and the corresponding optimal PMU locations 
can vary greatly. 

In this paper, we focus on the problem of using PMU 
measurements of voltage phase angles for real time detection 
of line outages over wide areas (including the ones that 
occur at places where no PMU measurements are available 
at the control center). We assume that network state esti- 
mation during normal conditions is available over relatively 
slow timescales, (e.g., using NERC System Data Exchange 
(SDX) iflOl ). which provides the pre-outage base case system 
parameters and states. In flD, and 15), it was demonstrated 
that many line outages can be detected with only a subset of 
the buses providing PMU measurements. However, a com- 
prehensive understanding of the outage detection performance 
given arbitrary constraints on the number of PMUs to use is 
left open. 

We address this open question of characterizing the optimal 
trade-off between the number of PMUs used and the outage 
detection performance. The central problem is finding the 
optimal locations at which to collect PMU measurements, 
which is an NP hard combinatorial optimization. We consider 
both non-adaptive and adaptive PMU location selections. 

We first define a set of voltage phase angle signatures associ- 
ated with the potential outages. Intuitively, the outage detection 
performance depends on the degrees of separation among the 
outage signatures. We employ the minimum distance among 
the signatures as an indicator of their degrees of separation. 
Accordingly, we introduce the max-min distance criterion for 
optimizing the PMU locations given any constraint on the 
number of PMUs. A simple greedy heuristic is developed that 
provides lower bounds on the global optimum. Next, we show 
that finding the optimal PMU locations can be equivalently 
formulated as an integer programming problem (IP), which 
allows a relaxation as a linear program (LP) that provides 
upper bounds on the global optimum. We then develop a 
branch and bound algorithm for PMU location selection: at 
each iteration, we strategically fix a bus to either provide or 
not provide PMU measurements, and compute new lower and 
upper bounds using the greedy heuristic and LP relaxation. 

To understand the optimal trade-off between the number 
of PMUs and the minimum distance among the outage sig- 
natures, we apply the proposed branch and bound algorithm 
on detecting all single line outages in IEEE 14, 24 and 30 
bus systems ifTTl . and show that the globally optimal PMU 
locations can be found in a small number of iterations. The 
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optimal trade-offs show that having roughly one third of the 
buses to provide PMU measurements is sufficient for achieving 
the same outage detection performance as with all the buses 
providing PMU measurements. 

II. Problem Formulation 

We consider a power transmission network with N buses 
and L transmission lines. We denote the set of all the buses 
by TV = {1,2,..., N}, and the set of all the lines by C = 
{1,2,..., L}. We study the problem of line outage detection, 
where an outage event corresponds to the loss of a subset of 
lines in C. Clearly, there are in total 2 L — 1 different possible 
outage events, which correspond to all the non-empty subsets 
of C. In practice, however, there are several reasons for us to 
consider a much smaller subset of the outage events: 

• In case a single line outage results in the overloading 
and overheating of another line, it often takes minutes 
after the initial single line outage for the overheated line 
to trip. Thus, within a few tens of seconds, it is very 
unlikely to see many line outages happen together due to 
their probabilistic independence. 

> Not all the outage events have the same level of impact 
on the grid. For example, some line outages do not lead 
to the overheating of other lines, while others do. The 
latter have a greater impact on the grid as they may cause 
cascading failures. 

• As we would like to detect and locate outages in real 
time with good detection performance, the limitations of 
computational power may forbid us to consider too many 
outage events. 

Based on the above considerations and the real world situation, 
a set of outage events of interest shall be selected, which we 
denote by £ = {Ei,E 2 , Ek}- For example, in (3), J4), £ 
contains all single and double line outages, respectively. 

Given an outage event set £, we augment it by the non- 
outage event (i.e., normal condition) Eq, and denote the 
augmented set by £ = {Eq, E\, . . . , Ek}- We investigate the 
following problem: 

• The grid starts in normal condition Eq. 

• Either the grid stays in normal condition E = Eq, or an 
arbitrary outage event E = Ek G £ occurs. 

> We make a detection decision E G £ . The detection is 
successful if E = E, otherwise not. 

We detect outages by observing the voltage phase angles 
at a subset of the buses, where the phasor measurements are 
provided by PMUs. A system diagram for this problem is 
depicted in Figure Q] As outages shall be detected as quickly 
as possible to prevent cascading effects, we aim at detecting an 
outage within about ten seconds after its occurrence. We note 
that the phase angles typically stabilize within a few seconds 
after a line outage occurs and we assume that the power 
injections of the network remain the same within these few 
seconds. This is a reasonable assumption as power injections 
change relatively slowly. We make use of the stabilized phase 
angles before and after an outage to make detection decisions. 




Fig. 1. System diagram for the problem of line outage detection using 
PMU measurements of voltage phase angles. Buses and transmission lines are 
represented by circles and connecting segments. The solid circles are the buses 
where PMU measurements are communicated to the outage detection module 
in the control center. The solid links represent real time communications of the 
measurements. The dashed SDX / state estimation module and links provide 
the slow timescale updates of the system parameters and states to the outage 
detection module. The loss of line X is an example of single line outage. 

A. Pre-computing the Outage Signatures 

For each event Ek G £, we denote the corresponding 
stabilized phase angles at all the buses by 8^ G K Arxl . We 
employ the DC power flow model to compute <?( fe ) as follows: 

P = B (fe) (fc) ,fc = 0, 1,...,K, (1) 

where P G R Ar><1 denotes the power injections, and B^ k > G 
R NxN is determined by the line status information after event 
Ei- occurs lfl2l . We note the following properties of H^: 

B {k) l = 0, (2) 
and we can arbitrarily add the same constant phase angle 
to every entry of 6, without having any change of P. 
• rank(S' fc ') = N — 1 if and only if the grid is connected. 
In the case where the grid contains islands, rank(i?' fe )) = 
N — C where C is the number of connected components 
(islands) of the grid. 
We note that if an outage creates islands in the grid, the 
balance of the power injections may not be satisfied by the 
same P, and (fTJ may not hold anymore. For example, in 
Figure Q] the loss of line a will create a single bus island 
A which may not be self-balanced (i.e., the pre-outage power 
injection at A may not be zero). We restrict our consideration 
to the cases where the post-outage grid remains connected, 
and leave the islanding cases as future work. 

In practice, noisy versions of the pre-outage and post-outage 
phase angle vectors are observed: 

e = (o) + z°, e = e {K) + z\ (3) 

where E k ,k G 0, 1, . . . , K is the actual outage (or non-outage 
if k = 0) that occurred, and z ,^ 1 are observation noise 
vectors that account for the errors in the measurement data 
and the system parameter data. In the following sections, we 
assume that the noises at all the buses are i.i.d.. We will 
generalize this assumption in Section IV-AI The task of line 
outage detection can then be formulated as the following 
hypothesis testing problem: 

From observing the pre-outage and post-outage phase angle 
vectors 0, 9 (01, identify which event Ek G £ has occurred. 

For this hypothesis testing problem, 6^°\ 6^, . . . , can 
be viewed as the signatures of the events Eq, E\, . . . , Ek- 
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To identify any event Ek S £, all the signatures {6^ k \k = 
0, 1, . . . , K} shall be collected before an outage occurs. 
From (|T), we pre-compute the signatures by 

0« = B ( - k '> + P,k = 0,l,...,K, (4) 

where _B( fc ) + is the pseudoinverse of B^. We assume that the 
knowledge of the power injections P and the normal condition 
£?(°) matrix are available from system-wide data sources (e.g., 
the NERC SDX |[Tol . or other state estimation mechanisms 
that operate over relatively slow timescales, cf. Figure Q]) 
Note that, with the knowledge of the open line indices for 
event Ek, the post-outage B^ matrices (fc = 1,2, . . . , K) 
can be derived from the normal condition JE?(°) matrix. As a 
result, while the grid is working under normal conditions, the 
signatures {9^ k \k = 0,1,..., A'} can be pre-computed by 
(|4]i in preparation for detecting the potential line outages in £. 
Intuitively, the more separated the signatures are from each 
other, the better detection performance can be achieved. In 
later sections, we characterize the separation by the minimum 
distance among the signature set. 

From |0, a phase angle vector remains functionally 
equivalent by adding the same constant to all its entries. 
Without loss of generality (WLOG), we can choose any one of 
the N buses as the reference bus, and subtract the phase angle 
at the reference bus from the phase angles at all the N buses, 
such that the phase angle at the reference bus is always kept 
at zero. In principle, it does not matter which bus among the 
N buses is chosen to be the reference bus. However, it does 
matter in the case when only a subset of the buses' phase 
angles are available. This is because once we choose a bus as 
the reference bus, we implicitly assume that the phase angle 
at this bus is available. 

B. PMU Location Selection 

We now establish the main problem of this paper, namely, 
PMU location selection. In (01, synchronous and timely up- 
dates of the complete phase angle vector require all the 
buses to have PMU measurements, and moreover low latency 
communication links that convey the PMU measurements from 
all the buses to the control center that performs timely outage 
detection. 

In practice, however, it may not be economically desirable 
that every bus has a PMU and constantly communicates the 
PMU measurements to the control center, particularly when 
there is very little performance loss with only a fraction of the 
buses providing PMU measurements. We motivate the selec- 
tion of a subset of the buses to provide PMU measurements 
in the following two application scenarios: 

1) As the cost of PMU installation in high voltage transmis- 
sion networks is relatively high, we may want to install 
PMUs only at a subset of the buses to reduce cost. 

2) Suppose over time the cost of PMUs drops and all (or 
many) buses have PMUs installed. Due to the infor- 
mation redundancy in the PMU measurements, it may 
not be effective to communicate the data from all the 
buses to the control center, consuming an unnecessary 
amount of the expensive low-latency communication 



link capacity, and also incurring a longer data collection 
delay. Thus, the control center may select only a subset 
of the PMUs to provide measurements. 
In the following sections, we address the second application 
scenario, in which the PMU location selection can adapt to the 
changes of the power injections P and the network B matrix 
over time. We show later in section IV-BI that the results can 
be directly extended to the first application scenario where a 
non-adaptive PMU location selection is required. 

The problem of PMU location selection entails the follow- 
ing two questions: 

• Given that we want to choose M (2 < M < N) of 
the N buses to provide PMU measurements, which M 
buses should we choose, and how can we characterize 
the corresponding outage detection performance? 

• What is the optimal tradeoff between the number of PMUs 
M and the outage detection performance? 

As mentioned in the last section, it makes no sense to have 
a bus chosen as the reference bus without having a PMU 
measuring its phase angle. Therefore, to find the optimal PMU 
locations, we cannot simply choose a reference bus arbitrarily, 
as it does lose generality. Instead, we traverse all the N cases 
of choosing each bus as the reference bus. 

From now on, we assume that we have chosen a reference 
bus r during this outer traversal, and optimize the location 
selection of the other M — 1 PMUs. We assume that all the 
phase angle vectors have been adjusted so that their r th entries 
always equal 0: 

6»< fc) =0,Vfc = 0,l,...,A. 

Accordingly, we denote by Ai the subset of buses with PMU 
measurements in addition to r, Ai C {1, 2, . . . , N}, \Ai\ = 
M—l. As a result, for any two phase angle vectors, we cannot 
distinguish their entries at the other N— M+l buses in Af\Ai. 
Therefore, for any phase angle vector 6, we project it into a 
sub-vector by extracting the M — 1 entries of 6 whose indices 
are in Ai (it is not necessarily an orthogonal projection) and 
denote the projected vector by 6m £ R( M_1 ) xl . Accordingly, 
every set Ai leads to a set of projected signatures {0 M \ k = 
0,1,..., A'}. 

Clearly, with different choices of Ai, the degrees of sepa- 
ration among the projected signatures can differ considerably. 
Intuitively we want to optimize Ai to get better separation 
among the projected signatures, and thus better detection 
performance. 

We model the separation among the projected signatures 
(and hence the outage detection performance) by the minimum 
distance in p-norm among them: 

d min (M) = mm - e M \, (5) 

0<z<.?<-K 

where p is a parameter to choose, and p = 2 corresponds to 
the Euclidian distance. Given M as the total number of PMUs 
to use, Ai can then be optimized under the following Max-Min 
Distance Criterion: 

max d min (Ai). (6) 

M,\M\=M-l,r£M 
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Clearly, (O is an NP hard combinatorial optimization of the 
set M. 



C. An Integer Programming Formulation 
Note that d m i n (M) can be re-written as 

d mia (M) = (mm(w T ®))> 

where 



w G 



oJVxl 



and Wj 



1, if i G .M or i = r, 
0, otherwise. 



(7) 



(8) 



• G R v 2 J, and its columns are constructed by 



collecting the following ( K ^ 1 ) N x 1 vectors 
V0 < i < j < K, 



(9) 



where the | • \ p operation is applied elementwise. 

In other words, the indices of the non-zero entries of w 
denote the M buses that are chosen to provide PMU mea- 
surements (including the reference bus r). We name the binary 
vector w the bus selection indicator vector. Consequently, © 
is equivalent to the following integer programming problem: 

maxmin(ii> T 0) 



s.t. Wi G {0, l},i= l,...,N, 
Wi = M, w r = 1. 



(10) 
(11) 



JV 
i=l 



Clearly, finding the global optimum of dTOb requires a worst- 
case computational complexity of 

III. PMU Location Selection with Max-Min 
Distance Criterion 

In this section, we provide two algorithms for solving the 
combinatorial optimization of PMU location selection ©. 

A. A Greedy Heuristic 

We develop a greedy heuristic that generates a series of 
PMU location selection solutions M2, M3, . . . , Mn for M = 
2,3, ... ,N respectively, that satisfy the following consistency 
property: 

M 2 C M 3 C . . . C M N - (12) 

We present the greedy algorithm as follows in a slightly more 
general form with an arbitrary initial set of selected buses 

Mini'. 

Algorithm 1: Greedy PMU Location Selection 

Given M and an initial set of buses Mini'- 

m = \M ini \,M m = Mim. 
Repeat 
m <— m + 1, 



M m = M m -i U < argmax d min (M m -i U {n}) 
Until m = M. 



(13) 



a starting set of chosen buses Mini, we choose another 
M — lA^ini] buses one by one: At each step, we keep the 
already chosen buses; from the remaining buses, we choose 
the one that maximizes the current step's minimum distance, 
and include it in the set of the chosen buses. When the only 
prior knowledge of the bus selection is the chosen reference 
bus r, Mini — {r}, and Algorithm 1 generates the set of 
greedy consistent solutions (fTSl i as M increases from 2 to N. 

B. A Branch and Bound Algorithm 

First, we note that the integer programming formulation (TlOb 
has a concave objective function, and hence has the following 
relaxation as a convex optimization: 

maxmin(u; T 0) (14) 



s.t. < Wi < 1, i 

N 



1. 



,N, 



(15) 



Wj 



M, w r = 1. 



In fact, this convex optimization can be equivalently cast as a 
linear program lfT3l . Accordingly, the optimal value of < TT~4b 
serves as an upper bound, denoted by U\, on the global 
optimum of (TlOb . Meanwhile, Algorithm 1 with Mini = {r} 
provides a lower bound, denoted by L\. 

Remark 1 (A note on the rounding heuristic): Another 
lower bounding heuristic is to find an integral solution by 
rounding the fractional solution obtained from the relaxed 
problem ( fl4l i. In particular, we consider the heuristic by 
rounding the M largest fractional entries to 1, and the others 
to 0. In all scenarios that we simulated, we compared this 
rounding heuristic with the greedy heuristic: somewhat to 
our surprise, the greedy heuristic uniformly outperforms 
the rounding heuristic. The reason is that, in the relaxed 
fractional solution, it is often some very small non-zero 
entries that are critical in the sense that losing them will 
drastically reduce the minimum distance. Consequently, for 
lower bounding the global optimum, we propose the greedy 
heuristic instead of the rounding heuristic, as the former 
is both much cheaper computationally and much better in 
performance. Thus, we use the relaxation technique not to 
provide a fractional solution to round, but to upper bound 
the global optimum and to develop a branch and bound 
algorithm that can significantly improve the greedy solutions 
in a few iterations. 

For any bus n ^ r, dTOb can be split into two sub-problems 
by fixing w n to be either or 1: 

max min ( it; T ) 



s.t. Wi G {0, 1}, i=l,. 

N 



,N, 



Wi = M, W n = 0, W r = 1. 



i=l 



and 



maxmin(u; T 0) 

w 

s.t. Wi G {0, 1}, i = l,. 



,N, 



(16) 
(17) 



(18) 
(19) 



N 



In other words, given the total number of PMUs M and 



^ Wi = M, w n = l,w r = 1. 



i=i 
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Similarly to ([141 ). relaxations of these two sub-problems can 
be formed by replacing ( fTTI i and ( fT9b with <JT5j, and they 
provide upper bounds, denoted by and on the global 
optimum of ( fT6] l and ( fl~8T > respectively. Meanwhile, applying 
the greedy heuristic under the constraint w n = or w n = 1 
provides lower bounds, denoted by ^ an d 4* > on these sub- 
problems' global optima. Define 

U 2 = max{4 0) , u ( 2 1] }, and L 2 = max{4° \ l { 2 1] }. (20) 

Then, U% and L2 are new upper and lower bounds on the 
original global optimum dTDl [14|. 

More generally, the above splitting procedure with relax- 
ations and greedy heuristics can be applied on the sub- 
problems themselves to form more children sub-problems with 
upper and lower bounds. For example, for any bus s, (s ^ 
n, r, ) ( [ToT l can be further split into two sub-problems by 
adding yet another constraint w s = or w s = 1 respectively. 

We define the following upper and lower bounding oracles, 
as well as an oracle that returns the next bus to split: 

Definition 1: Oracle UB(C) takes a constraint set C as in- 
put, where C specifies a set of buses whose selection indicator 
variables are pre-determined to be either or 1. An IP under 
the constraints C is formed, a relaxation is solved, and the 
optimum of this relaxation is output by UB(C) as an upper 
bound on the optimum of the constrained IP. 

For example, in ( TToT l and ( fT8l l. the constraint sets are 
C(°) = { Wn = 0,uv = 1} and C (1 > = {w n = l,w r = 1}, 
respectively. 

Definition 2: Oracle LB(C) takes a constraint set C as input. 
An IP under the constraints C is formed, a greedy solution is 
found by Algorithm 1, and the achieved objective value is 
output by LB(C) as a lower bound on the optimum of the 
constrained IP. 

Definition 3: Based on the order of the buses chosen by 
Algorithm 1, Oracle next(C) outputs the first bus that is 
chosen by this heuristic. 

When a sub-problem with constraints C needs to be split 
further, next(C) is the bus we choose to perform the splitting 
by fixing w next (c) to be either or 1. 

We now provide a branch and bound algorithm as in 
Algorithm 2 where i max is the maximum number of itera- 
tions allowed. As the algorithm progresses, a binary tree is 
developed where each node represents a constraint set. The 
leaf nodes are kept in S. The tree starts with a single node, 
{w r = 1}, corresponding to the prior knowledge that the 
reference bus is r who uses a PMU. When a sub-problem 
corresponding to a leaf node C* is split into two new sub- 
problems, the two new constraint sets and become 
the children of the parent constraint set C*. 

In Algorithm 2, (1221 is a generalization of d20l i. It means that 
the current global upper bound equals the highest upper bound 
among all the leaf node constraint sets. This is true because all 
the leaf nodes S represent a complete partition of the original 
parameter space fBl . At the beginning of every iteration, in 
choosing which leaf node to split OTT i. we select the one that 
gives the highest upper bound (i.e., the current global upper 
bound). It is a heuristic based on the reasoning that, by further 
splitting this critical leaf node, a lower global upper bound 



may be obtained, (whereas splitting any other node will leave 
the global upper bound unchanged.) At iteration i, the current 
upper and lower bounds on the global optimum are available as 
Ui and L{. When these two bounds meet, i.e., Ui — Li < e, the 
solution that achieves the current lower bound is guaranteed 
to be globally optimal. 

Algorithm 2: 
PMU Location Selection using Branch and Bound 

Initial step: i = 1, 

the initial constraint set: C% = {w r = 1}, 

the initial set of leaves of the tree of constraint sets 

(initially a single node): S = {Ci}. 

Compute Ui = UB(d), L x = LB{d). 
While Ui — Li > e or i < i max , repeat 

Choose which leaf node constraint set to split: 

C* = argmax{f7B(C)}. (21) 
ces 

Choose the next bus to split, n = next(C*), 
Form two new constraint sets, 

C{°\ = C* U {w n = 0}, Cgi = C* U {w n = 1}. 
In the set of leaves S, replace the parent constraint set 
C* with the two children c|°\ and C^\: 

S^(S\{C*})U{Cl Q + \}U{C^ + \}. 
Compute new upper and lower bounds for the two new 
constrained IP: 

UB{cfl), UB(CH\), LB{Cfll LB(Cll\), 
Update the global upper and lower bounds, 

U i+1 = max{[/B(C)}, L t+1 = max{£B(C)}. (22) 

i 4- i + 1. 

Choose the best achieved solution so far: 

C = argmax CeS LB{C). 

Return the greedy solution under the constraint set C. 



We note that Algorithm 1 is a degraded version of Algorithm 
2 with just one iteration. As the total number of possible 
constraint sets is 2^ (corresponding to the 2 N bus selection 
indicator vectors w), Algorithm 2 is guaranteed to converge 
in 2 N iterations (and in practice much less as will be shown 
later.) To limit the algorithm's run time, a maximum number 
of iterations i max can be enforced as in Algorithm 2. 

Finally, we define ^achieve to be the number of iterations 
used to achieve the globally optimal solution, and i pr ove the 
number of iterations used to prove its global optimality. In 
other words, it takes i ac hieve iterations for the lower bound to 
reach the global optimum, while it takes i proV e iterations for 
both the upper and lower bounds to reach the global optimum. 
As will be shown next, typically we have i ac hieve *C iprove- 

IV. Performance Evaluation and the Optimal 
M-d min Tradeoff 

In this section, we simulate the proposed algorithms in 
IEEE 14, 24, and 30 bus systems using the software toolbox 
MATPOWER |H). We set the outage event set £ to be 
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IEEE 14 bus system 




Number of PMUs 



Fig. 2. Tradeoff between the number of PMUs and the maximum achievable 
d m i n among the projected outage signatures; IEEE 14 bus system. 



IEEE RTS 24 bus system 




Number of PMUs 

Fig. 3. Tradeoff between the number of PMUs and the maximum achievable 
d m i n among the projected outage signatures; IEEE RTS 24 bus system. 

the set of all single line outages that do not create islands. 
Using Algorithm 2, we find the globally optimal solutions for 
PMU location selection, and characterize the optimal tradeoff 
between the number of PMUs used and the minimum distance 
in Euclidean norm among the projected outage signatures. 

We denote the maximum achievable minimum distance as 
a function of M by 

d* min (M),M = 2,3,..., N. (23) 

We plot d* min (M) in Figure El and ffl for the 14, 24 and 30 
bus systems respectively. In comparison, the greedy solutions 
achieved by Algorithm 1 are also plotted. 
We make the following observations: 

. Having M w f PMUs is roughly sufficient to achieve 
the same optimal d m - m as with M = N (i.e., all the buses 
having PMU measurements). 

• As M increases from 2, every additional PMU signif- 
icantly increases d* min {M), until roughly d^ in (N) is 
achieved (with M « y as mentioned above). 



IEEE 30 bus system 




Number of PMUs 

Fig. 4. Tradeoff between the number of PMUs and the maximum achievable 
d m i n among the projected outage signatures; IEEE 30 bus system. 

• When the available number of PMUs is sufficiently large 
(M greater than -j ~ in these examples), globally 
optimal performance can be achieved by simple greedy 
solutions. 

The number of iterations needed for Algorithm 2 to reach 
the globally optimal solutions, i a chieve and i proV e, are summa- 
rized in Table 1, 2 and 3, for the 14, 24 and 30 bus systems 
respectively. 

We make the following observations on the efficiency of 
Algorithm 2: 



Table 1: Number of iterations to reach the global optimum, 14 bus. 
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Table 2: 


Number of iterations to reach the global optimum, 24 bus. 
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Table 3: 


Number of iterations to reach the 


global optimum, 30 bus. 
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x -| q~ 3 Simulation results for Algorithm 2, IEEE RTS 24 bus system, M = 7 
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Fig. 5. The instant upper and lower bounds on the maximum achievable 
minimum distance as Algorithm 2 iterates; IEEE RTS 24 bus system, M = 7. 



In these three systems, Algorithm 2 always finds the 
globally optimal solution in less than 19 iterations. For 
solving the M = 4 case in the 30 bus system, i ac hieve = 
19. For most other cases, i ac hieve is much smaller. 
When M is relatively small, it takes a much larger 
number of iterations to prove that the solutions achieved 
within 19 iterations are indeed globally optimal. The 
maximum i pro ve is 395 for solving the M = 7 case in 
the 24 bus system. For this particular case, the upper 
and lower bounds achieved as Algorithm 2 progresses 
are plotted in Figure 

Achieve and ipro V e are much smaller than the combinatorial 
complexity of this NP hard problem. E.g., for solving 
the case of M = 7 in the 24 bus system, an exhaustive 
search has to traverse ( 2 7 ) = 346104 PMU location se- 
lections, whereas Algorithm 2 finds the optimal solution 
in ^achieve = 17 iterations, and verifies its optimality in 
*Di-ove = 395 iterations. 



The fact that i ac hieve is typically small demonstrates that 
Algorithm 2 is very efficient in finding the globally optimal 
PMU locations. We note that, in practice, the sometimes larger 
«prove does not matter at all. This is because we always use the 
best solution found within a predetermined run time, and as 
long as iachieve is sufficiently small, we will find the optimal 
solution (albeit without proving its optimality). 

Finally, as an example of the optimal PMU location se- 
lection, we depict the buses that are chosen for the 30 
bus system with M = 10 (cf. Figure [6j. In this case, 
bus 1, 5, 8, 9, 14, 21, 22, 24, 26, 29 are chosen to provide PMU 
measurements, achieving the same minimum distance among 
the outage signatures as with all the 30 buses chosen. As 
shown in Table 3, since i ac hieve = 1 for M = 10, this globally 
optimal solution can in fact be found by the simple greedy 
method of Algorithm 1. 




Fig. 6. One line diagram of the IEEE 30 bus system, and the optimal PMU 
location selection with M = 10. The buses enclosed by rectangles are the 
optimally chosen buses. 



V. Discussion 

A. Normalizing Distances by Noise Variances 

The max-min distance criterion we have used is for coun- 
tering the effect of the noise (fJJ by separating the outage 
signatures as much as possible. In general, the variances of the 
noise in the PMU measurements may differ at different buses. 
This shall be taken into account when defining the effective 
distances among the outage signatures. 

For two outage signatures and denote their dif- 
ference in absolute value by AO = |0W — 6^>\, where | • | 
is applied elementwise. The n th entry A9 n is the distance 
between the two signatures at the n th bus. 

Suppose now that the noises are independently Gaussian 
with variances al,a%, . . . ,ajj at the N buses. We can then 
normalize the distance at the n th bus by a factor of 



A0 = £ -1 A0, 



(24) 



where S = diag(<7i, C2, . . . , oat). After this normalization, it 
is now appropriate to add up all the buses' contributions in 
separating the signatures by computing the p-norm of A6: 



N 



IA0II 



(25) 



Accordingly, the definition of the minimum distance among 
the signatures (|5]l is generalized as follows: 



AM) — min 

0<i<j<K 



Q(i) _ gU) 



M 



8 



For the IP formulation ([Tol l, it was assumed that £ = al for 
some a in the previous definition of (O. With an arbitrary 
S, the matrix of the distances among the outage signatures is 
then generalized by the following normalization: 

= S p ®. (26) 

B. Extension to Non-Adaptive PMU Location Selection 

We have provided two application scenarios in Section Hl-BI 
that motivate the problem of PMU location selection. In the 
previous sections, we have focused on the second scenario 
in which the outage signatures k = 0, 1, ... , A}, and 

hence the optimal PMU location selection, depend on the 
current network setting P and B. 

In the first scenario, however, location selection for PMU 
installation is a planning problem whose solution must work 
non-adaptively in the sense that the installed PMUs are not 
transferable to other locations. As a result, the PMU locations 
must be chosen to accommodate, if not all, the most typical 
network settings in terms of P and B. As we consider wide- 
area transmission networks, this requirement is addressed as 
follows: 

• The network topology B is in general slowly and slightly 
changing, and an approximate estimate of it would suffice 
for computing typical sets of outage signatures. 

• Typical power injections P are available from grid statis- 
tics for high voltage transmission networks: collecting 
the typical data for different seasons and during daytime 
and nighttime would suffice for computing typical sets of 
outage signatures. 

With a set of typical power injections Px, P2, . . . , Pt, and 
an estimated B, there are T typical sets of outage signatures: 

{0 (M) , k = 0, 1, . . . , A'}, t = 1, 2, . . . , T. (27) 

The minimum distance (O can be generalized as 

d min (M)= min min - 0%%. (28) 

l<t<T0<i<j<K Jl Jl 

Similarly to ©, we construct © g R ArxT '( A 2 +1 ) i and its 
columns are constructed by collecting the following T ■ ( t 1 ) 
A x 1 vectors: VI < t < T, V0 < i < j < A, 

\g(.i,t) _ gU,t)\P } (29) 

where the | • \ p operation is applied elementwise. 

The problem of non-adaptive PMU location selection with 
(|28T > can then be formulated as the same integer programming 
problem ([Toi l, for which Algorithm 1 and 2 can be applied. 

VI. Concluding Remarks 

We have studied the problem of PMU location selection for 
detecting line outages in wide-area transmission networks. We 
have first introduced the voltage phase angle signatures of the 
outages. We have then established the PMU location selection 
problem as maximizing the minimum distance among the 
set of outage signatures, and have shown that it can be 
equivalently formulated as an integer programming problem. 
For this NP hard problem, we have developed a greedy 



heuristic and a linear programming relaxation, providing a 
lower and an upper bound on the global optimum respectively. 
Based on this heuristic and relaxation, we have proposed 
a branch and bound algorithm to find the globally optimal 
PMU locations. Using this algorithm, we have characterized 
the optimal trade-offs between the number of PMUs and the 
minimum distance among the outage signatures in IEEE 14, 
24 and 30 bus systems. For all the simulated cases, the optimal 
PMU locations are found in at most 19 (and in most cases a 
much less number of) iterations. From the optimal trade-offs, 
we have observed that it is sufficient to have roughly one third 
of the buses to provide PMU measurements in order to achieve 
the same outage detection performance as with all the buses 
providing PMU measurements. 

For future work, while we have employed the max-min 
distance criterion among the outage signatures, it is interesting 
to investigate the use of other criteria in characterizing outage 
detection performance. We would also like to examine the 
optimal PMU location selections with outage event sets more 
general than line outages. Finally, it remains an interesting 
open question how to take into account the needs from 
multiple tasks simultaneously (e.g., outage detection and state 
estimation) while selecting PMU locations. 
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